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A general nonperturbative theory of the low-energy electron propagator is developed and used to calculate 
the single-particle density of states in a variety of systems. This method involves the decoupling of the electron- 
electron interaction through a Hubbard-Stratonovich transformation, followed by a saddle-point approximation 
of the remaining functional integral. The final expression is found to be the tunneling analog of the infrared 
catastrophe that occurs in the x-ray edge problem; here, the host system responds to the potential produced 
by the abrupt addition of an electron during a tunneling event. This response can lead to a suppression in the 
tunneling density of states near the Fermi energy. This method is adaptable to lattice or continuum models of any 
dimensionality, with or without translational invariance. When applied, the exact density of states is obtained 
for the Tomonaga-Luttinger model, and the pseudogap of a fractional quantum Hall fluid is recovered. 
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I. INTRODUCTION 

In condensed matter, tunneling experiments provide a rich 
source of information. From the atomic scale analysis of a 
scanning tunneling microscope to the tunneling behavior of a 
macroscopic supercurrent, the tunneling behavior probes in- 
trinsic properties of the sample and provides opportunities to 
test theory against experiment, as well as generate new ques- 
tions. In the case of local tunneling, the tunneling current, or 
more specifically the conductance, is governed by the local 
electronic density of states (DOS) of the system of interest; 
thus, DOS measurements by tunneling experiments have be- 
come a common and indispensable tool of condensed matter 
physics. Therefore, it is of great importance to be able to the- 
oretically describe the tunneling behavior of such systems. To 
make this connection from the theoretical side, one has to cal- 
culate the relevant electron propagator of the system; as the 
DOS is related to a single-particle Green's function. 

Besides standard diagrammatic techniques, there are many 
analytical and numerical methods available to calculate the 
DOS, such as density functional methods,' which can reli- 
ably provide the DOS over a wide-energy range for weakly 
or non-correlated Fermi liquids. When coiTelations or 
electron-electron interactions become strong other techniques 
exist, such as bosonization- and density matrix renormal- 
ization group"* (DMRG) for one-dimesional systems and 
more recently dynamic mean field theory'' (DMFT) and its 
extensions"' in higher dimensions. Although, there have been 
tremendous strides in the area of strongly coiTelated elec- 
trons, there are physical systems where all current methods 
breakdown. Then one is usually left only with exact diago- 
nalization, which by its very nature is severely restricted by 
the exponentially increasing size of the Hilbert space. One 
such system is the edge of a quantum Hall bar,^^-^ *^ where the 
experiments of Grayson et al.,"^ on the tunneling spectra at 
the edge of a sharply confined two-dimensional quantum Hall 
system, are still theoretically unresolved. Thus, it is important 
to develop new methods and new insights into the tunneling 
behavior of such systems. In the follow sections we outline 
and demonstrate a novel approach to accurately calculate the 



low-energy DOS in a wide variety of cases, from Fermi and 
Luttinger liquids to the bulk fractional quantum Hall system. 
Ultimately it is hoped to apply the method developed here to 
the quantum Hall edge, where existing theoretical methods 
fail. 

In this article we derive an expression for a general 
fully interacting single-particle Green's function by perform- 
ing a Hubbard-Stratonovich transformation, decoupling the 
electron-electron interaction, followed by a saddle-point ap- 
proximation of the resulting functional integral. The final ex- 
pression provides a physically appealing picture that brings 
together semi-classical charge spreading theories'" and x- 
ray edge physics, "''"'''^''^ '^ as related to tunneling into low- 
dimensional and or strongly coiTelated systems. Although, an 
exact expression for the Green's function, within the saddle- 
point approximation, is obtained, an exact evaluation is gener- 
ally not possible. The main obstacle to this is determining the 
saddle-point itself, followed by evaluating the needed electron 
correlation function for the saddle-point field. Thus, we solve 
for an approximate saddle-point and then proceed to evaluate 
the coiTelation functions in the presence of this approximate 
field: either exactly or by cumulant resummation of a pertur- 
bation series. Within these approximations we find both qual- 
itative and quantitative agreement when compared to known 
results and other methods. It is also shown in Sees. Ill and IV 
that this approach absorbs our previous, physically motived 
but somewhat ad hoc, approaches''^''^''**''''' to this line of in- 
vestigation, by deriving them as limiting cases of the method 
presented here. Finally, in Sec. V we apply the full formalism 
to calculate the tunneling DOS of a bulk fractional quantum 
Hall system, where we recover the experimentally observed 
pseudogap in the DOS. 



n. SADDLE-POINT APPROXIMATION 

We start by considering a general D-dimensional interact- 
ing electron system, possibly in an external magnetic field; the 
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second quantized grand-canonical Hamiltonian is 
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where 11 = p+ ^ A, and vo{r) is any single-particle potential, 
which may include a periodic lattice potential or disorder or 
both. Here, e is the magnitude of the electron charge, i.e., 
e = |e|. Apart from an additive constant we can rewrite H as 
Ho + V, where 



and 



V : 



d"r' 6n{r) U{r - r') (5n(r'). 



(3) 



Hq is then the Hamiltonian in the Hartree approximation. The 
single-particle potential u(r) includes the Hartree interaction 
with a self-consistent density nQ{r), 



v{r)^v,{r)+ / dVC/(r-r')no(r'), 



with 



.„(r) = ^(vl/t(r)vl/,(r))^, 



(4) 



(5) 



where the chemical potential in Hq has been shifted by 
-C/(0)/2, and (O)o = Tr(e-'3^°0)/Tr(e-'^^°) denotes an 
expectation value with respect to the Hartree-level Hamilto- 
nian. In a translationally invariant system the equilibrium den- 
sity is unaffected by interactions, but in a disordered or inho- 
mogeneous system it will be necessary to distinguish between 
the approximate Hartree and the exact density distributions. 
The interaction in (3) is written in terms of density fluctua- 
tions; 



6n{r) = Y,'^Ur)^,ir) - no{r). 



(6) 



We want to calculate the time-ordered Euclidean propaga- 
tor (h = 1) 

G(rfaf,riCTi,To) - -(r«',,(rf,ro)^..(ri,0))^ (7) 

for large tq or equivalently at low energy. From which, the 
local DOS 

A^(r , = - - V Im G{r<T, rcr, iuj ^ uj + iO+) (8) 



is obtained by analytic continuation of the Fourier transform 
of the local Green's function. In the interaction representation 
with respect to Hq, 

. (Tvl/,^,(rf,ro)^<,,(ri,0)e-/o''<^-v'(r)^ 
G(rfcrf,ricri,To) = 



(9) 



Next we introduce a Hubbard-Stratonovich transformation of 

the form^^"* 



■ ^ / (5n U dn 



D(j> e 



(10) 



where the auxiliary fields are real-scalar functions, i.e., (j) : 
M^+^ M. To make the transformation well-defined the 
fields should satisfy 

j d^rdTd^r'dr' (P{r,T)U'\rT,rT')(j){r',T') > -oo. 

(11) 

Here, the inverse denotes a functional inverse; 
/ f{^, i)r\t y) = - y) and {/(r, r) = U{v)5{t). 
This Hubbard-Stratonovich transformation can be under- 
stood, from a quantum electrodynamics (QED) point of 
view, as a reintroduction of photons as the mediator of the 
electron-electron interaction in the c ^ oo limit. Using (10) 
in (9) gives 



G(rfcrf,riCri,ro) = TV 



where 
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(12) 



g(rfcrf,ricri,To|0) = 

- (r*^^. (rf , ro)^,, (ri, 0) ^ d-r^(r,r) Snir.r)^^ 

(13) 

is a noninteracting correlation function in the presence 
of a purely imaginary scalar potential i(j){r,T). M ~ 

(Texp(— dr V))q^ is a constant, independent of tq, and 
will remain an unknown overall prefactor. 

Because the bare interaction U in (1) is spin-independent, 
the spin dynamics remain trivial; therefore, spin labels will be 
suppressed, except where needed for completeness or clarity. 
The inclusion of spin related interactions within this formal- 
ism is possible. One simply needs to construct the appropriate 
spin-dependent Hubbard-Stratonovich transformation. 

Equation (12) can be written as 



G(rf,ri,ro) ^ Af 



D(f>e 
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Del) e 



(14) 



with an action 



^ J dTdT'd°rd^r'(j){r,T)U-^{r^r',T-T')(j}{r',T') 



lng(rf,ri,To|?!)). 



(15) 



Although, mathematically unjustified in the absence of a large 
or small parameter we evaluate the functional integral by a 
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saddle-point approximation. We define the saddle-point This requires the saddle-point field to satisfy 



as the field where the action is stationary, i.e., = 



</'sp(r,T) 



d^r' C/(r - r')('r«'(rf,To)^'(ri,0)5n(r',r)e 



g{'ri,ri,To\(j)sp) 

Expanding the right side in powers of (j>sp gives, 

0,p(r, r) = y d^r' C/(r - r') Ci(r', t) + 2J d'^T'dr" C2(r'r, r"r")(/<sp(r", r") + • • • 

where Ci and C2 are the first and second cumulants of a cumulant expansion of (13) defined by 

5(rf,ri,ro|0) = Go(rf,ri,To)exp j Ci(rr)(/.(r, r) + Cslrr, rV)0(r, T)0(r', t') + 

and are completely giving in terms of noninteracting Green's functions Gq{y, t); 

.Go(rf,r,To - r)Go(r,ri,r) 



Gi(r,T) 



Go(rf,ri,ro) 



and 



G2(rr,rV) = Jx(r, r', r, r') - lGi(r, r)Gi(r', r') - ^ 

2 2 2Go(,rf,ri,Toj 



(16) 



(17) 



(18) 



(19) 



Go(r, r', T - r')Go(r', Ti, T')Go(rf , r, tq - r) 
f Go(r', r, t' - r)Go(r, r;, r)Go(rf , r', tq - r') 



(20) 



r 



where x(r,r',T, t') = — {T6n{r,T)6n{r' ,t')) is the 
density-density correlation function. 

The exact solution of (17) can formally be written as 



(r, r)^tj d^'r'dT' U,s{rT, r'T')p(r', r'), (21) 



where 

[/eff(rT,r'T') 



[[/-i(r - r', T - r') - 2G2(rr, r'r') + •••]" 



is a highly non-linear inverse operator, and 

= -iGi(r,T). 



(22) 



(23) 



Expanding the action, Eq. (15), about the saddle-point field to 
second order and performing the functional integral"* gives 

G(rf,ri,ro) = ^(rf , r;, ro)5(rf , rj, To|(?!)sp), (24) 

as the saddle-point approximation for the interacting Green's 
function, where 



contains the fluctuation determinant. It is believed, but left un- 
proven, that this prefactor has weak tq dependance and there- 
fore will be neglected in any applied calculation. Equation 
(24) can also be written as 



G(rf,ri,ro) = ^(ro)Zsp(To)Gsp(rf , ri, ro)e' 



(26) 



where Zsp(ro) = (re'-/''^^'^''''*=p(""'^)"(''^^))^^, and the 
saddle-point Green's function 

Gsp(rf,ri,To) = 

\ / Ho 



(27) 



satisfies the Dyson equation 

G,p(r,r',T,T') = Go(r,r',T-r')- 

^ J dT"ciVGo(r,r",r-r")<A.p(r",r")G.p(r",r',r",r'). 

(28) 

Alternatively, using a cumulant expansion for 
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g(rf,ri,To|0sp), 
g(rt,ri,ro|0sp) = G'o(rf , ^i, tq) 

along with the formally exact expression for the saddle-point, 
Eq. (21), the interacting Green's function in the saddle-point 
approximation can also be written as 

G(rf,ri,ro) = ^(To)Go(rf , rj, tq) 

^ g— 1/2 / drdr'd^rd^r' p(r,r)Ueii{rT,r'r')p{r' ,t') (30) 

Equation (26) or equivalently (30) is the main result of this 
work and the starting point for the following sections. 

Although, formally exact, within the saddle-point approxi- 
mation, an exact evaluation of (26) or (30) is beyond reach, as 
it involves finding the solution to a non-linear integral equa- 
tion of infinite order for the saddle-point field, Eq. (17). Even 
so, making the simplest of approximations for Ucff, or equiv- 
alently (j>sp, captures a great deal of nontrivial physics. In the 
following sections we analyze both expressions [Eq. (26) and 
(30)] under various approximations and explore the underly- 
ing physics of the results. 

III. CUMULANT APPROXIMATION OF SADDLE-POINT 
EQUATIONS AND CHARGE SPREADING 
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FIG. I: (Color online) Normalized space-time contour plot of the 
charge density, Eq. (23), for a 1-D electron gas with rt = ri = 
and Tgep^ = 10. Notice, at the points (ri,0) and (rf,ro) it is 
highly localized, in accordance with the addition and removal of an 
electron. The stripes or bands are caused by Friedel-like oscillations. 
See Ref. [17] for further details. 



From a pure aesthetic point of view the cumulant represen- 
tation, Eq. (30), seems to have the most physically appeal- 
ing form. Within the suggestive notation of (30), it implies 
the effect of interactions is to renormilize the noninteracting 
Green's function by an action term e^^''^'^°\ where 



(31) 

that describes a space and time dependent charge density p 
interacting through f/cff. In fact, in Ref. [17] a very similar 
expression was found, although arrived at from an almost en- 
tirely different approach. There it was shown that p behaves 
similar to what one could loosely interpret as the charge den- 
sity of the tunneling electron, or more precisely the added and 
removed electron associated with the Green's function (7). In 
other words, p is highly localized around the space-time points 
(rj, 0) and (rt, tq) and in between it de-localizes, as it relaxes 
and interacts with the system. As an example. Fig. 1 shows 
a space-time contour plot of this "charge density" for a one- 
dimenstional electron gas. 

Returning to Eq. (30), as has been mentioned and can be 
seen, finding the exact effective interaction Ucs is unobtain- 
able. So one is force to make approximations. The simplest of 
which is replacing the effective interaction by the bare inter- 
action Ucs — > U, this neglects the higher cumulants and other 
terms of (22) [or (16)]. This actually involves two separate 
approximations for Eq. (30), one for the saddle-point and the 
other for g(rf , ri, To|0sp)- This is because to obtain (30) we 



assumed knowledge of the exact solution for the saddle-point 
field. Upon making an approximation for the saddle-point the 
term e""^"^ becomes a selective resummation of perturbation 
theory for (/(rf , ri, ro|(/)sp). Ultimately, this is what was ef- 
fectively done in Ref. [17] and the failures of the formalism 
presented there, specifically in the quantum Hall system, can 
be traced back to these two approximations. Nevertheless, it 
was also shown that even within these relatively simple ap- 
proximations the formalism was powerful enough to recover 
Fermi liquid theory in 2- and 3-D systems, while also obtain- 
ing the exact DOS of the 1-D Tomonaga-Luttinger model. 

If needed, as in Sec. V, to go beyond this two pronged ap- 
proximation one should solve for the full g{rf, r^, To|(/)sp), for 
an approximate (p^p- Doing this and tunneling's connection to 
the infrared catastrophe is developed next. 



IV. CONNECTION TO THE X-RAY EDGE PROBLEM AND 
THE X-RAY EDGE LIMIT 

Here, it is shown that the saddle-point field is a generalized 
x-ray edge potential. An x-ray edge potential is a localized 
time dependent field; generally of the form 

0,r(r,r) = [/(r)e(To-T)e(T). (32) 

An infrared catastrophe can be caused by the singular screen- 
ing response of a system to potentials of the form of (32). This 
behavior is known to be responsible for the singular x-ray op- 
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tical and photoemission spectra of metals,"'-"'' Anderson's 
orthogonality catastrophe, '^^''^ and even the Kondo effect. 

The x-ray edge problem, the singular absorption spectra of 
metals, was predicted by Mahan," while studying "excitons" 
in metals; the hole of the electron-hole pair is created by the 
absorption of a soft x-ray. It is assumed the hole has a large 
effective mass and thus the Fermi-sea of the metal responds to 
a potential of the form of where the hole is generated at 
T by the absorption of a photon and then later recombines 
at r = tq. 

Note, if — i0xri apart from phase factors, Eq. (26) is 
equivalent to the x-ray edge problem,-" where Ggp contains 
the "exitonic" part and Z^p is the Anderson orthogonality con- 
tribution. This connection between tunneling and the x-ray 
edge field can most easily be seen by starting with the approx- 
imate (for transparency only) saddle-point field 



d°r'U{y-v')p{v',T). 



(33) 



As was covered in the previous section, p plays the role of an 
effective charge density of the tunneling electron. Expanding 
it in multipole moments^^*" 

p(r,T) =p"(T)(5(r)+pi(T)r-V(5(r) + -- - , (34) 

where p"(t) — J dPr r"p{r, t). It can easily be shown'^ that 
p" = Q{tq ~ t)Q{t). By keeping only the monopole term, 
the saddle-point field is then 



0,p(r,r) «i[/(r,0)e(ro-T)e(r). 



(35) 



This is the x-ray edge potential, which corresponds to the po- 
tential produced if the tunneling particle had an infinite mass. 
Indeed the saddle-point field is a generalized (p^^, that ac- 
counts for both dynamic interactions and recoil of the particle. 

We define the so-called x-ray edge limit as the limit where 
the approximate saddle-point field is taking to be of the form 
of (35). Thus, the fully interacting Green's function in this 
limit is 



G(rf,ri,To) 5(rf,ri,To|i0xr)- 



or 



G(rf,ri,To) ~ Zsp(ro)G'sp(rf,ri,To)e 



/ 4>xrno 



(36) 



(37) 



In this limit, apart from energy shifts, tunneling exactly re- 
duces to the x-ray edge problem. This x-ray edge limit was 
studied in detail in Refs. [16,18,19] for a 1-D metal and the 2- 
D quantum Hall system, where qualitative results were found. 
Of course it is desirable to go beyond this limit and include 
the full dynamics of the saddle-point (within a given approxi- 
mation of the integral equation for ^sp)- This can be difficult 
even in the simplest systems; one commonly has to find nu- 
merical solutions to the many integral equations, if possible. 



V. APPLICATION TO A FRACTIONAL QUANTUM HALL 
FLUID 

Because the kinetic energy, and hence the relaxation and 
recoil of a newly added electron, in the lowest Landau level 
(LLL) of a quantum Hall system is completely quenched, one 



would expect this system to be the prototypical example of an 
infrared catastrophe occurring during tunneling. In the LLL it 
is known experimentally and theoretically-^'^"*'^^'^^'^^'^^'-^'-'"''" 
a pseudogap develops in the DOS near the Fermi energy. In 
Ref [16], it was shown in the x-ray edge limit, 0sp 0xr, a 
pseudogap is indeed recovered, but this limit was introduced 
by hand, based only on plausible physical arguments. It is not 
obvious using the full formalism presented here that such a 
result can be obtained. 

Here, we apply the method of the previous sections to a 
two-dimensional electron gas in a quantizing magnetic field. 
Choosing B = —Be^, and the symmetric gauge A = 
{By/2)ex — {Bx/2)ey, in the zero-temperature limit the non- 
interacting Green's function for the LLL is 



Go(r,r',T) = r(r,r')[j^-e(T)], 



(38) 



where 



9.TV ^ 



{rrj 



27r ^ 2™m! 

m=0 



.gim(e-e')g-rV4g-r'V4 



(39) 



and 0„j(r) are the single-particle eigenfunctions (in units 
where the magnetic length I — yJhcjeB = 1), < < 1 is 
the filling factor, and 0(t) is the Heaviside step function. The 
charge density, Eq. (23), for this system is therefore 



27r 



-[e(To-r)e(T)-H. 



(40) 



In systems such as this, with ground-state degeneracy, 
the leading order approximation for (/igp, (/)sp(r, r) w 
i \ dPr'XJiv — r')p(r'.r), fails the integrable condition, 
Eq. (11), for any meaningful interaction in the /3 — > oo limit. 
Therefore, one has to find another approximate solution of 
(17) by including higher-order terms, in </)sp, to obtain a suit- 
able saddle-point. At a minimum one has to solve the lin- 
earized version of (17). This is what we know turn to. Addi- 
tionally, we will assume a finite inverse-temperature /3, only 
taking /3 — > cx) at the end if desired. At finite temperature the 
LLL Green's function is 

Go(r,r',T) = r(r,r')K(eo - p) - e(r)]e-(^o-^)^ (41) 

where riyiuj) is the Fermi distribution, eQ is the single-particle 
energy of the lowest Landau level, and the chemical potential 
p is defined such that 



C30 

E 



1 



expL 



m)] + i 



(42) 



Note, the Green's function is aperiodic, Go{T—f3) = —Gq{t), 
as it should. The finite-temperature charge density is then 



Hr|V2 



2tt 



-[e(ro-T)e(r)-nF(eo-M)]- (43) 



From (17), the linearized integral equation for 0sp, setting 
Tf = Ti = 0, is 
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,p(r,r)=i j d\'U{v-v')p{y',T) + J d\'Sr" j dr' U {v - v')x{y' ,v" ,T,T')cj,,^{v" ,t') 



+ 



d^r'd^r" / dr 



.(eo - m) - 0(r - T')]K(eo - ^) - e(ro - r)] - K(6o - a^) - e(r' - r)]K(eo - - e(ro - r')] }</>.p(r", r'), 

(44) 



r 



where x(r,r', r,T') = — " (2^Tr)" ^ In the low- 

temperature limit, the solution is (see appendix A) 

<^,p(r,r) ^ ^ 1 dV'[/(r-r')e"l"-'l'/2e(ro-r)e(r). 

(45) 

Because this is a valid saddle-point field, we can now deter- 
mine the interacting Green's function. Within the cumulant 
approximation. Sec. Ill, using (45) gives 



G(0,0,ro) 



i/- 1 

( 

2tt 



(46) 



where A = (47r)-i // dVd2r'[/(r - rOe^l'-'I'/^e-klVz. 
Equation (46) leads to a delta function or "hard-gap" in the 
DOS 



iV(0, uj) = const X S[uj - A(l - v)], 



(47) 



instead of the observed pseudogap. Mathematically a pesudo- 
gap can be modeled as a function that decays faster than any 
power-law at low energy. To obtain this behavior implies that 
in the time domain the function should exhibit a similar de- 
cay rate as r — > cx). Therefore, the fully interacting Green's 
function would be expected to be a nonalgebraic and or possi- 
bly nonanalytic function of tq. The cumulant approximation, 
which is basically a selective resummation of perturbation the- 
ory in (j)sp or more explicitly for this case perturbation in Atq, 
isn't powerful enough to capture such behavior; an exact solu- 
tion of the Dyson equation for the saddle-point Green's func- 
tion, Eq. (28), is required for this system. 



With (45), the Dyson equation for the saddle-point Green's 
function is 



G,p(r,r',r,r') = T{r,r')[np{eo - - e(r - r' 



'^lp-(eo-A')(^— r') 



'0 

^ J dTidVid2r2r(r,ri)[nF(eo-Ai)-e(r-ri)]e-(^o-'^)(^-^i)f/(ri-r2)e-l"-^l'/2Gep(r2,r',Ti,T'). (48) 



The solution of which is 
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G,p(r, r', r, r') = ^ 0„(r)0;^(r')a„(T, r'), 



(49) 



where 

/ (np(eo - Ai) - l)B(r - r') + np(eo - /x)e(T' - r)exp(-A„To) 



exp[-A„(r - r')] exp[-(eo - Ai)(r - r')] 

(50) 



and 



^ — j d^rd^r' 4>* {v)U{v ~ v')e^^^'^^ /"^(j) (r) (51) tion in the LLL. For a screened Coulomb potential of the form 



are the diagonal matrix elements of the saddle-point interac- 



U{y) = e 



2exp(-a|r|) 



(52) 
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FIG. 2: The DOS (up to an overall constant prefactor) of the LLL. 
Found by analytic continuation of (55), using the experimentally 
relevent parameters /? = 10, a = .05, and v = 1/3, in units of 
eV^and L 

oo 

X^ = e' fdk'^^^t^^F,il + m,l;-ky2l (53) 
J va'^ + 



where ^F^{a,b;z) is the confluent hypergeometric function 
of the first kind. Experimentally the screening length a would 
typically be the distance to the nearest gate or to the second 
2-D electron gas in bilayer systems. 
From appendix B 

oo 
m=0 

(54) 

This finally gives the interacting Green's function as 

G(0, 0, To) - Z,p(ro)G,p(0, 0, to) exp (i/a" Vq) 

^ OO 

= 27r n ~ ^F^'^o - m) + ^¥(^0 - m) exp(-A™ro)] 

m— 1 

X {npie^ - A*) - l)cxp(-AoTo)exp (t/a^Vo). (55) 

Figure 2 shows the DOS obtained from (55), by analytically 
continuation using the maximum entropy method. As can be 



seen, a pseudogap is present in agreement with experiment. 
Although, the width of the peak is of the order seen experi- 
mentally, the actual location cannot be totally trusted. This is 
because we have solved for Ggp using a long time approxima- 
tion for Gq. Then used the solution to find Z^p, (see appendix 
B) but short time information enters into Zsp. This is a well- 
known result from the x-ray edge problem, where unknown 
energy shifts can occur. The inclusion of higher Landau lev- 
els in Go would fix this. 

VI. CONCLUSIONS 

By starting with an exact functional integral representation 
of an interacting Green's function, followed by a saddle-point 
approximation, we have shown that the low-energy tunneling 
behavior of electronic systems is governed by the response of 
the host system to the time-dependent potential produced by 
the newly added electron. In the x-ray edge limit the infrared 
catastrophe occuiTing during a tunneling event can be exactly 
mapped to the well-known x-ray edge problem. Within vari- 
ous approximations for the saddle-point we have obtain quali- 
tative as well as quantitative results in a wide spectrum of sys- 
tems including: Fermi liquids, the Tomonaga-Luttinger liq- 
uid, and quantum Hall system. 

Future work will involve applying and comparing this 
method to other exactly solvable models, such as the 1-D Hub- 
bard model and ultimately to the edge of a quantum Hall fluid, 
where existing theoretical methods fail. 
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APPENDIX A: SADDLE-POINT FIELD FOR THE 
QUANTUM HALL SYSTEM 

Here, we solve for the saddle-point field including the lin- 
ear term in (f)sp of the integral equation (17). This is needed 
because the simplest approximation for the saddle-point ob- 
tained by keeping only the zeroth order term of Eq. (17) is not 
in the integration domain of the functional integral. For the 
quantum Hall system this then reads 
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</),p(r, T)=ij dPr' U{r - r')p{r', r) + j dPr'dPr" j dr' U{r - r')x(r', r", r, rO<Asp(r", r') 



+ J^J d^'cPr" j dr'[/(r-r')e-l-'lVVI-"l^/2{[e(To-r)e(r)-nF(eo-M)][e(To-r')-nF(eo-/x)] 



- K(eo - m) - e(T - r')]K(eo - m) - e(ro - r)] - K(eo - m) - e(T' - T)]K(eo - m) - e(To - r')] }</.sp(r", r'). 

(Al) 

Using x(r,r',T,T') = g-l—'lV^ and 

e(To - r) [e(To - r') - e(r - r')] - e(T' - r)e(To - r') = (A2) 

gives 

K{r,r)=i J d\' U{r - r')p{r\ r) - ' ^^{2^)2"^^'° ~ / ^^VdV'^ dr' C/(r - r')e-l-'-"l^/Vsp(r", r') 



^ nAe^-^l)[l^~Me,-^i)] j ^2^,^2^n j dr' U{v - v')e-\^'\'^e-\^"\l^4>.^{v" ,t'). (A3) 



Introducing the Fourier transforms 

<^,p(k, zn„) = Jd^rJ dr </.sp(r, r)e'^''-e-''^-- and <A,p(r, r) = J (Cf^ '^^p^^' ^^")e"'''""e''^-', 


(A4) 

where if2„ is a bonsonic frequency, the integral equation in Fourier space is 

</.,p(k,zf]„) = zUik)p{k,zn„) - %(^o - /^)[;^-^f(^o - m)] ^(k)e-l'-lV^^,p(k,^»„0/3'^^a„„o'^»o„.o 

+ ~ ^'^1',;)^^'° ~ ""^^ ^(k)e-l-l^/^ / rf^fc'e-l'''l^/^0.p(k',ifi„O/35.a„.o<5.a„,o. (A5) 

Thus for iQn 7^ 

0sp(k, iQr, ^ 0) = iC/(k)p(k, iQ„ ^ 0) (A6) 

and for if2„ = 

</.,p(k,iQ„ = 0) = it/(k)p(k,in„ = 0) - j£^n^{^o - - ^F(eo - /x)]C/(k)e-l'^l'/203p(k,^^^„ = 0) 



/3 



(eo - - nF(eo - ii)]U{\<i)e-\''\"/^ j rf2A;'e-l'''l'/Vsp(k',in„ = 0). (A7) 



(27r) 
or 

0.p(k,iQ„ = 0) = , m^^)pO^.^^^ = ^) 

' l+(Jpnp(eo-/z)[l-np(eo-/x)]C/(k)e-|k|V2 
/? nF(eo-/x)[l- 71^(60 -M)]t/(k)e-l'-lV^ 

(27^) = 



+ ^ n (eo -/.)[! -n,(eo-M)]t/(k)e-l--l/2 . rf2^,^-|k'|V2<^ = 0) (A8) 
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This is a Fredholm integral equation with a degenerate kernel whose solution is simply given by^^- 

' 1 + (2^«F(eo - - «F(eo - A*)]C/(k)e-|kP/2 

+ _0) P np(.o^M)[l~^F(eo-M)]t/(k)e-l'-lV^ 

^ " " ^ (2-)^ 1 + - - n,(eo - M)]f/(k)e-l>^IV2 ^^'^^ 

where 

i[/(k)p(k, 0)e-l'^l'/2 



A(iU.n = 0) = — , „ , ■ (AlO) 

P f n^{e,-^,)[l-n^{e,-^i)]U{\.)e-\^\'/^e-\^\'l^ 



Converting back to r space 

0sp(k, ^) - ^ E «f/(k)p(k, zf]„)e-'""^ - /3-iiC/(k)p(k, if]„ = 0) + /S" Vsp(k, = 0), (Al 1) 

where the iO, — Q term has been subtracted from the sum over frequencies. In the /3 — > oo limit 

0,p(k, r) ^ -Lf/(k)e-l''l'/2e(To - r)e(T), (A12) 

or in real space 

,(r, r) = ^ y f/(r - Y')e'\^'\^ '^Q{to - t)Q{t). (A13) 

I 



This is the form of the saddle-point field that is used in Sec. V. constant integration trick found in Ref. [12]. 

From the linked cluster theorem, 

APPENDIX B: SADDLE-POINT S-MATRIX r ^ ^ , . / > 

Here, the saddle-point scattering matrix Zsp is found from 
Gsp. This is a straightforward generalization of the coupling where 



= I] T / dTld'^ri ■ ■ ■ dnd^n 0sp(ri, n) • • • 0sp(r;, r,) (Tn(ri, n) • • • n(r,, ri))differentconnected (B2) 



M = 7 / dTid'^ri ■ --dTid^ri (/)sp(ri,Ti) • • • (?!)sp(r,, Ti)G'o(ri, r;, ti, r,) • • • Go(r,, ri, r;, Ti+). 



(B3) 



Changing the summation limit leads to 



°° 1 /■ 

M T— . / dTid^ri ■ ■ ■ dr;+id^r,+i0sp(ri,ri) • • • 0sp(r;+i, r,+i)Go(ri, r/+i, ti, t,+i) • • • Go(r;+i, ri, t;+i, r+). 

1=0 

(B4) 

Introducing 

G«p(r, r', T, T+) = Go(r, r', T,T+)-t^J dr"d^r" Go(r, r", r, r")4p(r", r")G«p(r", r', r" , r+) (B5) 



and using Eq. (B4) can be written simply as 

1 

TTi 



^ = d^^' (B6) M^i j dT j d'^rj^ d^(/.,p(r,T)G«p(r,r,T,T+). (B7) 
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Reinstating spin for completeness gives 



M^iY,JdrJd''rJ^ 0,p(r, r)G«p(r, r, r, r+, a). 



(B8) 
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/ 



where A is a symmetric matrix. 
^* In general ijisp is a complex valued function, : R^+^ C, 
and is no longer in the domain of integration. Technically one has 
to deform the functional contour to pass through the saddle-point 



hy (j) ^ + "^sp and then expand the integrand around = 0. 
SeeRef. [19] for details. 

This can only be done in the absence of ground-state degeneracy. 
In such cases one has to first solve Eq. (17) and then expand (/>sp 
in multipole moments. For example, see Sec. V. 



